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High-sensitivity 27 A1 nuclear magnetic resonance (NMR) measurements of aluminum metal under 
high pressures of up to 10.1 GPa reveal an unexpected negative curvature in the pressure-dependence 
of the electronic density of states measured through shift and relaxation, which violates free electron 
behavior. A careful analysis of the Fermiology of aluminum shows that pressure induces an electronic 
topological transition (Lifshitz transition) that is responsible for the measured change in the density 
of states. The experiments also reveal a sudden increase in the NMR linewidth above 4.2 GPa from 
quadrupole interaction, which is not in agreement with the metal's cubic symmetry. 



PACS numbers: 76.60.-k, 62.50.-p, 71.20.Gj, 71.30.+h 

The fundamental properties of metals are well known 
PP, but nuclear magnetic resonance (NMR) at Giga- 
Pascal pressures supported by band structure calcula- 
tions open a new window for their deeper understand- 
ing. This will be demonstrated here with a surprising re- 
sult for an almost free electron system, simple aluminum 
metal - the most abundant metal of the earth's crust. 

Already in 1936, a decade before the discovery of NMR, 
Heitler and Teller [2] recognized that the relaxation rate 
(1/Ti) at which nuclei exchange energy with the lattice 
is much faster in metals than in non-metals due to the 
high density of electronic states (DOS) at the Fermi level, 
N(E F ). Later in 1949, Knight found a large NMR fre- 
quency shift (Knight shift) caused by electron spin para- 
magnetism [3] [4] , and Korringa [5] related 1/Ti to the 
temperature (T) independent Knight shift Ks- His fa- 
mous relation, T{FK 2 S = (j e /ln) 2 ft/(47r/c#), depends 
only on the gyromagnetic ratios of nucleus (7 n ) and 
electron (7 e ), apart from fundamental constants. Mate- 
rial properties may enter as electronic correlation effects 
through a correction factor [6] . 

Qualitative theories frequently describe N(E) with 
simple analytic functions. However, in real systems N(E) 
is always non-analytic at those energies where the band 
structure, T7(k), has extrema or saddle points. These de- 
tails can profoundly influence N(Ep), which is responsi- 
ble for many characteristic properties of a metal, includ- 
ing conductivity, specific heat, or the critical temperature 
of superconductivity. 

In an interesting scenario van Hove singularities in 
N(E) can be shifted by varying parameters like chem- 
ical composition, magnetic field, or pressure. If a sin- 
gularity passes through the Fermi level the metal un- 
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dergoes an electronic topological transition (ETT), of- 
ten referred to as Lifshitz transition, with a change in 
topology of the Fermi surface (FS) [7]. Unfortunately, 
changes in the chemical composition can induce disorder 
that masks subtle effects [8 , and the required magnetic 
field strengths often cannot be reached with normal mag- 
nets [9 j. Therefore, varying pressure appears to be the 
most appropriate tool for such studies (e.g., [TO]). 

The compressibility of Al metal [1] tells us that for a 
10% volume reduction, for which one might expect sub- 
stantial effects, we need a pressure of about 7.6 GPa. 
Such pressures can routinely be achieved with anvil cells. 
However, experimental methods that investigate the elec- 
tronic structure are rare under such conditions since in 
anvil cells two conical anvils with sub-millimeter culet 
size push on the tiny sample and the metallic gasket that 
encloses it. For example, in NMR experiments one needs 
to record a very weak, precessing nuclear magnetization 
with a resonant radio frequency (RF) coil that should 
fit the sample tightly, which was not possible with anvil 
cells, so far [TT] . 

A variety of different NMR setups have been employed 
in the past to investigate materials under such high pres- 
sure conditions [12-16]. Recently, a new high-sensitivity 
NMR anvil cell was introduced by some of us that uses 
an RF microcoil inside the high pressure region [17] , sim- 
ilar to designs used for low-frequency measurement of the 
susceptibility 18] and the de Haas- van Alphen effect [19] 
at high pressure. While this method has great poten- 
tial in other areas (as was shown, e.g., with the closing of 
the NMR pseudogap of a high temperature superconduc- 
tor [20]), its combination with bandstructure calculations 
seems to be particularly fruitful for the investigation of 
metals. 

Piston-cylinder type Moissanite anvil cells made of 
non-magnetic high-strength beryllium copper (BeCu) 
were loaded in Cambridge and readied for NMR in 
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FIG. 1: (Color online) 27 Al NMR of aluminum metal, (a) 
Spectra at 17.6 T as a function of pressure, (b) Pressure 
dependence of the total linewidth (squares) and the dipolar 
contribution (circles), (c) Field dependence of the FID at 
10.1 GPa (symbols) and ambient pressure (dashed line). The 
difference is defined as log[J(17.6 T)]-log[J(11.7 T)]. (d) FID 
and Hahn echoes of varying pulse separations (60, 80, 100 \is) 
observed at 5.4 GPa. 

Leipzig. Metallic aluminum powder (Alfa Aesar, purity 
99.97%, particle size less than 44 |xm) was filled into a 
9 to 10 turn RF microcoil wound from 12 |Ltm diameter 
copper wire that was placed inside the hole of the BeCu- 
gasket. The microcoil also contained a small ruby chip for 
pressure monitoring [21], Daphne 7373 and Glycerin were 
the pressure transmitting fluids. Various cells were used 
with home-built NMR probes that fit our standard wide- 
bore NMR magnets. Typical NMR experiments with tt/2 
pulse lengths of less than 1 ills were employed to excite 
and record free induction decays for shift and relaxation 
measurements. All measurements were carried out at 
room temperature. 

Some typical 27 Al NMR spectra up to 10.1 GPa to- 
gether with other experimental results are shown in 
Fig. [I] In panel (a) we see spectra as a function of pres- 
sure at 17.6 T. One observes a gradual shift of the reso- 
nance to lower frequency and an increase in the linewidth 
(and loss in intensity) as the pressure increases. The 27 A1 
nuclei in the metal are in a cubic environment and one 
expects only a single resonance line for this 7 = 5/2 
quadrupolar nucleus. At ambient pressure, we find a 
Knight shift Ks — (vo — Vret) / v ve { of 1640 ppm (the ref- 
erence was AICI3 • 71H2O), in good agreement with what 
has been reported before [22] . We find the free induction 
decay (FID) at ambient pressure to be nearly Gaussian, 
cf. dashed lines in Fig. [I] (c), with a second moment of 
9.86 G 2 , as expected for homonuclear dipolar coupling 
and in agreement with the literature value of 9.8 G 2 
[23] . (Beats in the FID at longer times, which were ob- 
served early on [23] and that are similar to those observed 
in CaF2 [24] can be explained in terms of properties of 
chaotic many-body quantum systems [25} [26]). 



In terms of the Korringa relation, we find at 295 K that 
T\T^ 1.85 sK, also in agreement with earlier studies that 
showed that the relation holds with the same constant 

I. 85 ± 0.05 sK between about 1 K and 1000 K [23]. 

In Fig. [I] (b) we show the total linewidth (squares) 
and the dipolar contribution 2a^/2 log(2)/7r to it (cir- 
cles), where a is the second moment determined from 
the beginning of the FID. As the pressure increases we 
observe that the dipolar linewidth does not change (from 
the change in lattice constants [27] we estimate a 10% 
increase), while the total width shows a rather steep in- 
crease above about 4 GPa to almost twice its ambient 
pressure value at 10.1 GPa. 

While no structural phase transition is expected to 
occur in our pressure range [28], we cannot exclude 
quadrupolar effects for 27 Al NMR, e.g., from mechani- 
cal strain. In fact, in such a case an increase in linewidth 
is expected. Furthermore, a shift of the resonance fre- 
quency might occur due to the quadrupole interaction 
impairing the exact determination of the Knight shift. 

In order to prove that the observed pressure-dependent 
shift is of magnetic origin we performed NMR experi- 
ments at various magnetic field strengths (B = 7.0 T, 

II. 7 T, and 17.6 T). Sample FIDs at 10.1 GPa are shown 
in Fig. [l] (c). We find that the shifts are indeed purely 
magnetic, cf. also Fig. [2] and that both linewidths con- 
tributions are independent of field. With the nuclear 
magnetic dipole broadening being independent of field 
we conclude that the additional broadening must be 
quadrupolar. We estimate the quadrupolar broadening 
to be about 13.5 kHz at 10.1 GPa, cf. Fig. [I] (b). In 
such a case one also expects that one can excite a weak 
Hahn spin echo that is absent for pure dipolar coupling 
[29] . Indeed, we find no Hahn echo at ambient pressure, 
but the appearance of a weak echo at high pressures, cf. 
Fig. [l](d). While we are not sure about the origin of 
this weak quadrupolar coupling at high pressures, now, 
we are certain that we can neglect it while interpreting 
the changes in shift and relaxation. 

The pressure dependence of the Knight shift, Ks(p), 
(also for different magnetic field strengths) is shown in 
the main panel of Fig. [2] Also shown (dotted line), is the 
expected dependence (Ks(p) oc N(Ep) ex V 2 ^ 3 ) for a free 
electron metal of volume V [27], which is clearly in con- 
tradiction with our measurements. Note that we do not 
expect the hyperfine coupling constant to change consid- 
erably with pressure as discussed below. The spin lattice 
relaxation at all pressures was found to be mono expo- 
nential, but pressure dependent. In the inset of Fig. [2] 
we plot the Korringa relation determined from our shift 
and relaxation data. It proofs that our observed changes 
in K and T\ can indeed be understood in terms of the 
hyperfine coupling to the s-electrons. 

In order to shed light on the discovered discrepancy be- 
tween free electron behavior of the shift and the experi- 
mental data, we decided to perform numerical bandstruc- 
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FIG. 2: (Color online) Experimentally observed Knight shift 
Ks at 11.7 T (triangles) and 17.6 T (squares). Also shown is 
the pressure dependence of Ks for a free electron gas (dashed 
line) and that predicted from the LDA calculations. (Inset) 
Pressure dependence of the Korringa factor l/(T\TK s ). 
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FIG. 3: (Color online) Electronic DOS of face centered cu- 
bic aluminum in LDA. (Main panel) Zero-pressure DOS, the 
Fermi energy Ef is indicated by a vertical line. (Upper in- 
set) Zoom near Ef for different pressures (see text). (Lower 
inset) Pressure dependence of Ns s (Ef), arrows indicate the 
positions of two ETT. 



ture calculations. The calculations were carried out with 
the full-potential local-orbital scheme, FPLO-9.00-34 [30] 
using the local density approximation (LDA Perdew- 
Wang 92) and a scalar-relativistic mode. The DOS was 
evaluated by linear-tetrahedron integrations with Blochl 
corrections on an extremely fine mesh of 300 x 300 x 300 
points in the full Brillouin zone (BZ) in order to achieve 
the required accuracy of the DOS close to singularities. 

Fig. [3] shows in the main panel the calculated total 
electronic DOS. From the band bottom to about -5 eV 
we find almost perfect free-electron behavior, but a num- 



ber of van-Hove singularities appear near the Fermi level. 
The pressure evolution of the total DOS close to E F 
with higher energy resolution is depicted in the upper in- 
set (the equation-of-state was obtained from a standard 
Birch- Murnaghan fit resulting in a zero-pressure LDA 
volume of 15.831 A 3 , a bulk modulus of 83.7 GPa, and its 
pressure derivative 4.4). In the occupied part two singu- 
larities move upward in energy with increasing pressure, 
while a weak singularity in the unoccupied part moves 
downward (arrows). In the lower inset, the 3s DOS at 
Ef (Ns 8 {Ef)) is shown as a function of pressure. Note 
that A^ 3s (£ 1 f) is non-analytic at those pressures where the 
topology of the FS changes and related van-Hove singu- 
larities pass Ey- The two visible critical points at about 
7.5 GPa and at about 25 GPa indicated by arrows are 
related to the two stronger singularities discussed above, 
while the third and weaker one is not resolved in this fig- 
ure (it almost merges with the singularity at 7.5 GPa). 

In order to compare the calculated DOS with the room- 
temperature NMR experiments, Ns s (Ep) was folded with 
a 300 K Fermi-Dirac function and the zero pressure point 
was normalized to Ks(0) (Fig. [2]). As one can see, the ex- 
perimental and theoretical data match almost perfectly. 
In the FPLO code, distortion of the wave functions un- 
der pressure is taken into account by including atomic- 
like wave functions beyond a minimum basis set. In the 
present case, the valence basis comprised 2s, 2p, 3s, 3p, 
3d, 4s, and 4p states. If a 4s contribution is taken into 
account, the high-pressure side of the theoretical curve 
shifts upward by about 6 ppm (in units of Ks). This 
means, we find an increase of the hyperfine coupling (av- 
eraged over 3s and 4s contributions) of 0.4% at 10 GPa, 
compared to its zero-pressure value. Semicore 2s and core 
Is contributions are found to be negligible. 

Qualitatively, the critical point at 7.5 GPa readily ex- 
plains the measured negative curvature of Ks(p): pres- 
sure drives the system through an ETT and the pressure 
dependence of Ns s (Ef) is stronger on the high-pressure 
side of the transition than on its low-pressure side. 

Fig. [4] illustrates the origin of the DOS singularities in 
the pressure dependence of the third-zone FS of Al. At 
p = the textbook four-armed "monster" is to be seen 
stretching along the square around X. It is characterized 
by three extremal orbits a, /?, and 7. At about 7 GPa, the 
monster has developed spikes which now touch the point 
W (arrow in Fig. [4] a saddle point of the band struc- 
ture at W crosses the Fermi level). Since another spike 
touches W from the neighboring BZ, the Fermi surface 
becomes multiply connected at this ETT [31] . While this 
topology change certainly has a drastic effect on mag- 
netotransport, it creates only a minor singularity in the 
DOS (arrows in the upper left insert to Fig. [3| and, thus, 
is not visible in the Knight shift. 

Only slightly above 7 GPa the monster breaks into 
parts. At this transition, the a-orbit vanishes and 
Ns s (Ef) shows a downward-kink (right lower insert to 
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FIG. 4: (Color online) Calculated Fermi surface of alu- 
minum. Blue (red) color indicates the occupied (unoccupied) 
side. Left part: The well-known third-zone electron monster 
at zero pressure. Right part: Detail of the FS with pres- 
sure dependence (indicated at the right-hand side). At zero 
pressure (lowermost figure), the detail is identical with the 
framed part of the left figure. Extremal orbits in the common 
notation [35] are indicated. The a-orbit vanishes at an ETT 
between 7 and 8 GPa and the 7-orbit vanishes between 24 
and 26 GPa. 

Fig. [3]). From this kink, Ks(p) obtains the observed 
negative curvature. Continued increase of the pressure 
lets the arms (or, spindles) around U shrink further un- 
til they vanish in another ETT. This happens at about 
25 GPa and causes an upward-kink in Ns s (Ef), which 
should manifest itself in a positive curvature of Ks(p) 
at p « 25 GPa that could not be reached with NMR 
experiments yet. 

Early studies of the FS of Al suggested that it is 
fairly well described by the free-electron model, see 
Ref. [32] and references therein. Later, de Haas-van 
Alphen experiments at pressures up to 0.7 GPa hinted 
at deviations from that model [33]. No such deviations 
were found in Fermi momenta measured by positron 
annihilation up to 10 GPa [34 , but the error bars 
associated with the reported data are sizable. Existing 
experimental information was used by Joss and Monnier 
to fit pseudopotential parameters and, subsequently, to 
predict the general linear strain dependence of Al FS 
areas [35] • In accordance with this parameterization, an 
ETT under uniaxial stress close to 0.5 GPa has been 
identified by magneto-thermopower measurements on 
Al whiskers [36 . Further, experimental data on Al 
alloyed with 6-10% Si have been interpreted in terms 
of an ETT [37]. We are not aware, however, of any 
previous experimental observation of ETT in Al under 
hydrostatic pressure. Our calculated transition pressures 
for the vanishing of the a- and 7-orbits are by factors 
of two to four lower than what one would estimate from 
the zero-pressure derivatives of the FS areas obtained 
by Joss and Monnier [35 . Thus, linear approximations 
for the pressure dependence of the FS area are obviously 
not in general valid for pressures beyond 1 GPa. 

To conclude, we have studied the Knight shift and re- 
laxation of aluminum metal as a function of pressure up 



to 10.1 GPa, extending the pressure range previously in- 
vestigated by a factor of 14 [38]. (Note that the authors 
in [38] found an increase in Knight-shift with pressure 
that we do not observe. However our first pressure point 
at 1 GPa is already higher in pressure than any of the 
data of the previous study.) Up to 10.1 GPa, we found 
a decrease of the Knight shift by 11% which is almost 
three times of what is expected for a free electron sys- 
tem. Since we also observed an increase in linewidths 
with pressure (above 4 GPa) that is inconsistent with 
homonuclear dipolar interactions we performed an ex- 
tensive study of the field dependence of the linewidths in 
order to exclude quadrupolar effects on the shift. We find 
the additional linewidth to be independent of field and 
we ascribe it to a small quadrupole interaction (13 kHz) 
that appears rather suddenly above 4 GPa. Since we used 
different anvil cells and different pressure transmitting 
fluids we believe that the measured, weak quadrupole 
coupling is not simply caused by strain, but a manifesta- 
tion of a charge density variation that breaks the cubic 
symmetry, as suggested earlier [39]; further experiments 
are currently under way to study this effect. 

Since the pressure dependent shift is magnetic and the 
Korringa relation does not change with pressure, we at- 
tribute the changes in Ks (p) and T\ (p) to changes in the 
electronic DOS. We performed numerical bandstructure 
calculations with high resolution and find quantitative 
agreement with the experimental data. We also find that 
the behavior is due to a kink in the electronic DOS at 
an electronic topological transition. While such transi- 
tions are normally only observable at low temperatures, 
our NMR data taken on a large enough pressure interval 
allow the detection of ETT even at room temperature. 
We predict that a further increase in pressure (up to 30 
GPa) will result in a further drop of the Knight shift by 
more than a factor of two. Finally, we expect a positive 
curvature of 27 K(p) around p = 25 GPa due to the dis- 
appearance of the spindle-like Fermi surface sheet close 
to the U-point. 
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